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ABSTRACT 

An experimental investigation of thermocapillary flow in droplets of a vegetable oil (partially 
hydrogenated soybean oil) immersed in silicone oil was conducted in a test cell with a heated top wall and a 
cooled bottom wall. The liquids are nearly immiscible and have equal densities at a temperature below the 
room temperature, thus providing a simulation of low-gravity conditions by reducing the buoyancy forces. 

The interfacial tension between the two oils was measured in the temperature range 20-50° C using a capil- 
lary tube and — — was determined to be negative. Droplets ranging in sizes from 3mm to 1cm diameter were 
dT 

injected into the silicone oil. The vertical temperature profile in the bulk liquid (silicone oil) produces 
temperature variations along the interface which induce variations in the interfacial tension. The flow inside 
the droplet driven by the resulting interfacial shear stresses was observed using a laser light- sheet flow 

visualization technique. The flow direction is consistent with the sign of — . The observed maximum surface 

velocities are compared to the theoretical predictions of Young et al. (1959). 


•NASA Resident Research Associate at Lewis Research Center. 


For short times after injection, the droplets were driven by this flow towards the hot wall above the 
matched-density temperatures until the droplet reached a point where the forces due to the flow and buoyan- 
cy were equal. After longer times, the droplets moved to the cold side due to suspected density changes 
caused by mass transfer from the droplets to the silicone oil. This was confirmed by tests under isothermal 
conditions, where it was observed that droplets of all sizes fell to the cold bottom eventually. 

Thus, even though the thermocapillary flow inside the droplets persisted for long times in spite of the 
mass transfer, the migration of droplets towards the hot side was eventually affected by uncontrolled buoyan- 
cy forces resulting from density changes due to mass transfer. While additional liquids are being tried, it is 

\ ■ > 

suggested from the present experience that reduced gravity experiments will probably be necessary to pro- 
vide unambiguous data for the migration of droplets. 

1 Introduction 

The investigation of the migration of droplets in liquids due to thermocapillary induced flow is sig- 
nificant not only in enhancing the understanding of such flows but also has direct application to material 
sciences and chemical engineering. It has been suggested that the low gravity environment available in space 
orbit can be exploited to provide advancements in materials processing (Neilson and Weinberg, 1977, Wein- 
berg, 1978). For example, containerless manufacturing of high technology ultrapure materials may be pos- 
sible in such an environment because of the elimination of heterogeneous reaction and nucleation effects at 
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the container walls, allowing larger undercooling of the melts, and substantially reducing the oftentimes 


deleterious effects of time-dependent buoyancy induced flows. However, such space manufacturing en- 
deavors involve complex, non-equilibrium systems oftentimes with multiple interfacial regions providing the 
potential for significant surface tension driven flows which could also be detrimental to product quality. 

Thus, it is imperative to develop an understanding of such effects through a coordinated earth-space re- 
search effort. The ground-based experimental study reported here was conducted to support such an effort 
in one special area — thermocapillary migration of immiscible droplets/bubbles in a continuous phase liquid. 

The migration of bubbles, or droplets, is an important aspect of many engineering sciences and space 
manufacturing issues. It is a prime candidate for many studies as it is without complications such as tree, or 
moving contact lines. In material science processes as in glass manufacturing, etc., gas bubbles can be 
formed from the byproducts of chemical reactions or gas trapped in the interstices of the raw material. In a 
1-g environment the migration of the bubbles due to buoyancy effects can be utilized as a separation 
mechanism. However, in a microgravity environment of space other forces must be utilized for the separa- 
tion. Young et ai (1959) pointed out that a fluid droplet placed in a thermal gradient environment will ex- 
perience a force due to the induced capillary flow. The interfacial tension will typically be minimum at the 
warm pole of the droplet, and will increase toward the cooler pole. The resulting tangential stress will drag 
the neighboring continuous phase fluid toward the cooler pole, and the reaction on the droplet will therefore 
be in the opposite direction. If the droplet is free to move, it will migrate toward warmer regions as a consc- 
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quence, and hence provide a separation mechanism sensitive to the magnitude and direction of the bulk 
temperature gradient. After the pioneering work of Young et al. (1959), the Marangoni convection of bub- 
bles in liquids was subject to many experimental investigations (Papazian and Wilcox, 1978, Hardy, 1979, 
Thompson, 1979, Nahle et al . , 1987). The investigations of the thermocapillary migration of droplets in liq- 
uid are more recent and less have been performed. Chifu et al. (1983) observed the flow on the surface of a 
droplet (freely suspended in a liquid) by injecting a small amount of a dyed surfactant into the interface. The 
velocity of the spreading surfactant front was measured and compared with the values resulting from their 
theoretical model. By also, using a two component system with a miscibility gap (methanol/cyclohexane) and 
cooling it down below the temperature where cyclohexane separates, they observed a thermocapillary migra- 
tion of the separated droplets in the direction of the warmer region. The thermocapillary migration of water 
droplets in butyl benzoate was investigated by Delitzsch et al. (1984). The two components of the liquid sys- 
tem that they used were immiscible and have the same density at about 32° C. Water droplets injected into a 
stably stratified layer of butyl benzoate sorted out according to their diameters in a manner that smaller 
droplets migrated to higher temperatures than the larger ones. The final droplet positions was determined 
by the balance of the thermocapillary (proportional to the surface area) and buoyancy forces (proportional 
to the volume). Seikmann et al. (1986) also used water droplets in silicone oil in a microgravity environment 
(D1 shuttle mission), but did not observe thermocapillary migration. Barton and Subramanian (1989), ob- 
served thermocapillary convection of ethyl salicylate droplets in diethylene glycol for Marangoni and 
Reynolds numbers much less than one. 
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In order to continue to increase the understanding of these phenomena it was decided that an important 


step would be to elucidate the flow inside and on the interface of a droplet, and its migration as well as its 
equilibrium position in a vertical temperature gradient. The visualization of the flow pattern would be an 
added confirmation that the mechanism for the migration is a result of the surface induced thermocapillary 
flow (in the opposite direction to the migration) and not due to any bulk induced motion (like buoyancy in- 
duced flow, for example, which would be in the same direction). Since conflicting results have been obtained 
using water systems, it was decided to use an easily available matched-density (at 15 C) liquid-liquid system 
consisting of vegetable oil as the droplet liquid and 5 CS. silicone oil as the continuous phase liquid. 

2 Theoretical Considerations 

To our knowledge, there is no theoretical treatment of thermocapillary migration of droplets for non- 
negligible Marangoni numbers. One choice for the purpose of comparison of experimental results is to refer 
to the classical linear theory of Young et al. (1959) as a guideline. They obtained the migration velocity V of 
the droplets as 


v P 2 £ (>-/>' )Q+/Q da_ Da 

6ft(2fi+fi') dz dT ( 2fi 4- 2>n ' ) ( 2a+ a ' ) 


( 1 ) 
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where D is the droplet diameter, g the gravitational acceleration, — the temperature gradient at infinity, 

the thermal interfacial tension coefficient, /taadft' the dynamic viscosities and a and a' the thermal dif- 
dT 

fusivities of the continuous phase fluid and the droplet liquid respectively. In the above equation and m what 
follows, a prime (’) designates the property in the droplet phase. The migration velocity is the sum of two 
terms, the first representing the velocity of a droplet due to buoyancy forces and the second term describing 
the velocity induced by Marangoni forces. When the droplet reaches a stationary position, Marangom and 
gravity forces balance each other and V = 0. 


k k’ 

Substituting/? = Pm [ 1 — ^ ( ^eq ~ ) ]> p ( ^cq T m )},a p'Cp 

(where k,C ? , k\ C p ' are the heat conduction coefficient and specific heat for continuous phase and the 
droplet liquids respectively) into the above equation, one can write. 


r„-7W(±) f 


( 2 ) 


where K = 




P m iU “-/*') ( 2+ ^") 


— is a material constant. Therefore the temperature dif- 


ference from the matched density temperature (Teq - T m ) at which the droplet is stationary (T eq ) is a 


6 


1 dT dT 

function of two controllable variables, namely — , and If is held constant, all droplets of equal size 

dT 

should experience the same deflection from the matched density temperature level. On other hand, if —j~ 

were to be increased, the above deflection would be larger. In the present study, experiments were performed 

dT 

to measure (T cq - T m ) withD and —as controlled variables. 

The maximum velocity on the interface according to YoungS theory occurs at the equator of the droplet, 
and is expressed as follows: 


Vmax. 


d o 

<-df> 


dT 

dz 


<§> 


0+/O (2+j) 


(3) 


This has been used for a comparison to our measured interfacial velocities on the droplet equator. 


3 Experimental Apparatus, Test Liquids and Procedure 


3.1 Experimental Apparatus 


The apparatus used in these experiments is shown in Figure 1. It consists of a rectangular cell filled 
with the test fluid with a movable top wall. The top and bottom walls are made from two copper plates, and 
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the side walls are made from transparent 1/2" thick plexiglass. Plexiglass has low conductivity and permits 
good flow observation. Copper tubes, of 3/8" diameter, shaped in a square coil are soldered outside the cell 
to both the bottom and the adjustable upper plates (milled and polished surfaces of l/4 ,! x5"x5", l/4 n x4"x4" cop- 
per plates) for cooling and heating purposes. The tubes of each plate are connected to separate constant 
temperature baths so as to adjust and maintain the temperatures of the plates. The four vertical walls are 
fixed to the bottom wall and sealed to prevent any possible leakage. The upper wall is attached from the top 
side on an adjustable column to enable changes in the distance between the copper plates. The upper wall 
also has two 1/4" diameter holes, one in the center and another 1" away from the center and side wall. These 
holes are used to fill the test cell with the continuous phase fluid, and also to provide access to position ther- 
mocouples for temperature measurements. The lower copper plate has a 1/8" diameter hole at the center 
with a drain valve attached to it from the outside. A compartment is, therefore, created between the copper 
plates. It is within this chamber that the experiments occur. In order to achieve a stably stratified liquid 
layer without any unwanted flows due to buoyancy, the upper plate was maintained at a higher temperature 
than the lower plate in all of the experiments. 


An injector was designed to produce droplets from 0.5 to 10mm diameter. The injector is a 0.25mm 
diameter and 3" long hypodermic needle and was connected, via a 1/8" diameter tygon tubing, to a syringe, 
which in turn, was supported by a holder at about the injection level. The syringe was manually actuated to 
produce the droplets. The injector hole in the side wall is placed at about the location of the matched-dcn- 
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sity temperature of the liquid-liquid system used in this experiment, and is sealed with a rubber cork. After a 
single droplet of a desirable size is generated the injector was very carefully pulled out of the side wall of the 
test cell to deploy the droplet. A laser light-sheet flow visualization technique has been used to study the 
flow field in the droplets. The optical components for this technique consist of a 0.5mm laser light beam 
passing through a cylindrical lens, thereby generating a sheet of light 0.3mm in width. The light-sheet can be 
used to illuminate flow tracer particles introduced in the flow field of interest. The motion of the particles (1 
micron aluminum oxide) reflecting the laser light can be observed and recorded on film and video cassette 
tapes for further analysis. This technique was used to observe the flow patterns and to estimate the droplet in- 
terface maximum velocities. 

3.2 Test Liquids 

Quite by accident we found that a nearly immiscible, non-toxic liquid which had a density very close to 
the 5 CS, silicone oil is vegetable oil. Therefore we chose this as the droplet phase. We measured most of the 
necessary liquid properties in our laboratory. The matched density temperature for the experiments reported 
here was measured in two ways: 1) Using a Mohr-Westphal specific gravity meter along with a HAAKE F4-C 
temperature bath, densities of both liquids vs. temperature were measured separately, and the best fit lines 
over the measured temperature range were obtained. The intersection of the two fitted lines was chosen as 

o c 

the matched-density temperature. This was at T m = 15 C where p m = 0.911 (also 
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5 = 0.001197 — -J- and 8' = 0.000703 — Jr were obtained fromthese measuremenls). 2) An 
r cm C cm C 

isothermal condition was created in two double wall containers, connected in series, by the same temperature 
bath. With the continuous phase liquid in the first container and the droplet liquid in the second, the bath 
was run until a predetermined isothermal condition in the entire volume of the liquids was achieved. The 
droplets were then placed into the continuous phase liquid. At temperatures not close to 

7m ( [ T - 7’ m | > or = 2° C ), the droplets either fell down to the bottom ( T> T m + 2° C ) or floated on top 

( T< T m - 2 °C ). For temperatures close to T m ( 1 7- T m \ < 1 °C, with T m = 15 °C), the droplets were 
neutrally buoyant, thus confirming the measurement of T m by the previous method. However, continuous 
mass transfer from the droplets to the silicone oil which occurred shortly after the droplets were inserted 
masked the precise determination of T m . In fact, in the vicinity of T m , peculiar results were observed 
wherein smaller droplets fell down to the bottom faster than bigger droplets, in contradistinction to the set- 
tling velocity of Hadamard and Rybczynski. Most likely, these results can be attributed to the mass transfer 
process. Mass transfer from the droplet oil to the silicone oil was also confirmed by observing precipitation 
in refrigerated samples of the silicone oil that came in contact with the vegetable oil. 

It was also noticed that the matched density temperature ( T m ) was changed by the aging of the 
vegetable oil over a few months, but the final experiments from which data is reported were conducted within 
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a five day period which, to our knowledge, did not seem to change the properties of the droplet liquid. 

T m = 15 °C was obtained at the beginning of this five day period. 

The viscosity of these liquids were measured in our laboratory at the room temperature by a Cannon- 

Fenske Routine viscometer. The interfacial surface tension was measured over a temperature range of 20 C 

to 50°C, by the use of a capillary tube (Rashidnia and Balasubramaniam, in progress). The data were suitab- 
ly fitted, and the thermal interfacial tension coefficient in the best fit line over the measured temperature 

range was calculated to be — 0.0088 An average thermal conductivity for the vegetable oil 

ri i 

( k' = 0.0003266 k ~ ) was taken from the Chemical Engineering handbook (Perry and Chilton, 

\ cm sec Cl 

1973). Thermal conductivity of the silicone oil (k = 0.00026 — ^ o r ) is given by DOW CORNING litera- 

\ cm sec e / 

ture. The Prandtl number and the thermal diffusion coefficient for the continuous phase fluid are 71.75 and 
2 

0.0007 . Table 1 summarizes the major properties of the fluid system used for this study. 

sec 

3.3 Procedure 

Before each experiment, the old liquid inside the cell was drained, and the cell was wiped, cleaned, and 
rinsed with fresh silicone oil. The cell (10x10x3 cm 3 test chamber) was then assembled and filled with 5 CS. 
silicone oil. Then, the circulating baths were turned on to bring the upper and lower boundaries of the cell 
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to the desired temperatures. This causes a vertical temperature gradient to be established in the cell. After 


waiting a sufficient amount of time (at least four hours) for steady state conditions to be achieved, the 
droplets were injected into the cell and their size, motion, and equilibrium positions observed through a 
Panasonic high resolution video camera and recorded on a super VHS video cassette. The diameter and posi- 
tion of the center of the droplets were also monitored by a cathetometer, and were manually recorded. 

The temperature in the host fluid was measured by a K-type thermocouple fixed to a 0.25mm moveable 
stem attached to a micrometer traversing mechanism. The thermocouple signals are read using a Keithly 871 
Digital thermometer. Data collected from experiments using three temperature gradients are reported in 
this paper. For each experiment, the period from starting the temperature baths to the end of data collec- 
tion lasted for about 9 continuous hours. The temperature gradient was checked at the end of each day's ex- 
periment and no measurable change (thermometer resolution = ±0.1 C) was detected. The temperature 
gradient was very linear in and beyond the range where the droplets migrated, except very close to the upper 
wall (within 1mm) where a small deviation occurred which could not affect the region of interest. The ther- 
mocouple was also traversed in two locations —in the center of the cell and half way to the side wall (2.5 cm 
apart). The temperature profilejin the cell, for the two locations, were very similar confirming a one-dimen- 
sional temperature field in the cell. 
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Figure 1 shows a schematic of the experiment and the optical arrangement utilized in the laser light- 
sheet flow visualization part of the experiment. Figure 2 is a photograph of the apparatus. The laser light- 
sheet flow visualization part of the experiment was conducted after the data collection was completed. 

Using the same technique, we were able to obtain average maximum velocities of the particles at the in 
terface of the droplets. The surface velocity has a maximum value at the equator (Young et «/., 1959) . 


In all the experiments described here, the droplet temperature at equilibrium is assumed to be the same 
as the temperature measured on an imaginary horizontal line passing through the droplet center (using the 
horizontal part of the cross hair of the cathetometer) that is one inch from the vertical axis of the cell (where 
the droplets were normally located). 

In order to obtain a measure of the magnitude of free convection flows in the cell, separate experiments 
were conducted. One micron size aluminum oxide particles, as tracers, were introduced into the continuous 
phase liquid. Sufficient time was allowed for the cell to return back to the expected steady state temperature 
gradient. Using a 10 mW power He-Ne laser light-sheet enabled visualization of particle motion in different 
vertical planes of the cell. Analysis of the video recorded images and long time exposure photography 

showed that free convection velocities were less than 10^ in the worst location (up to 1mm below the 

upper copper plate) which was along the surface direction (showing a thin cell far away from the center of 
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cell where the droplets were injected and reached their equilibrium positions). This minor velocity of tracer 


particles could be due to a small tilt angle of the upper wall and did not seem to affect the lateral position of 
the droplets that were observed for more than an hour after their injection. The rest of the test cell did not 

show any velocities of the tracers more than 5 which were mainly downward and that was attributed to 

sec 

gravity induced settling. 

When droplets were deployed, they were, for the most part of the experiment, located 3 to 10 radii from 
the horizontal walls, and 10 to 33 radii from the vertical walls of the test cell. The influence of the walls on 
the flow and the droplet equilibrium positions are discussed later. 

4 Results and Discussion 

All droplets analyzed were spherical, and the measured diameters and the equilibrium temperatures 
were constant during the time period of the experiment wherein the measurements have been made. A typi- 

D 2 

cal temperature versus normalized time, r = after droplet deployment is shown schematically in Figure 

3. t is the characteristic time for the droplets to reach their thermal equilibrium. After about 1.25tol.5r, 
droplets reached their maximum positions from the matched-density temperature T m . Observing this loca- 
tion for a short time (about 2— 3r) and making sure that it did not change, this location was taken as the 
droplet equilibrium position, and the necessary data were recorded. This process is shown schematically in 
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Figure 3. A slight plume was often observed behind the droplets resulting from dissolution of droplet liquid 
and the consequent changes in refractive index. Because of dissolution of some unknown component(s) that 
could change the density of the droplet liquid, the droplets and their central temperatures were tracked very 
carefully, and recorded up to the peak position shown in Figure 3. After the peak (longer times than 4r), 
droplets started to fall toward the cold wall. Typically, it took about 10r for the droplets to fall back to the 
T m level. 

Figure 4 displays the measured equilibrium temperature difference from T m versus the inverse of the 
droplet diameter for the three temperature gradients considered. Figure 5 shows the same data plotted as 

(-) 

( 7eq— T m ) versus - ^ - along with the predictions from linear theory (i.e. equation (2)). As can be seen 

from Figure 5, the measured data qualitatively follows the trend predicted by the linear theory, but is quan- 
titatively quite different. The droplets seem to have migrated farther away from Tm when compared to the 
predictions. Also, the data in Figure 5 does not collapse on a single straight line for all temperature 
gradients and droplet sizes, as anticipated by the linear theory. Table 2 shows the measured maximum inter- 
face velocities from particle tracing visualized by the laser light sheet and compared to the predictions 
from equation (3). The measured velocities are typically smaller by a factor of 2 to 3. Interestingly, this is 
similar to the measurements of migration velocities reported by Barton and Subramanian (1989). The flow 
Reynolds numbers for the present experiments were very low (i.e., Re< < 1). The Peclet numbers ranged 
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from 1 to 9. Figure 6 shows a long time exposure photograph of the largest section of a vegetable oil droplet 
showing the pathlines of the thermocapillary flow traced by the seed particles. 

While the droplets migrate farther away than predicted by linear theory (Figure 5), the interface 
velocities seem to be smaller (Table 2). Hence, stronger thermocapillary flow does not seem to be the reason 
for the excessive migration. By elimination, the larger migration distance must be because the effects of 
buoyancy are reduced compared to what linear theory anticipates. We speculate that this may be due to the 
following reasons: (i) The temperature distribution in and around the droplet (and consequently the density 
distribution and the buoyancy force) is different from linear theory predictions due to convection of energy. 
As mentioned before, the Peclet numbers range from 1 to 9. Effects due to energy convection are not negli- 
gible and the isotherms would be distorted compared to those in pure conductive heat transfer. Since we 
have not measured the temperature distribution in the presence of the droplets, we do not know what the 
isotherms look like and cannot evaluate whether buoyancy effects have increased or decreased compared to 
what linear theory anticipates. Also, to our knowledge there is no theoretical analysis predicting the 
temperature distribution for conditions of the current experiments, (ii) as mentioned before, the droplets at 
equilibrium positions were located 3 to 10 radii from the bottom wall of the test cell, with the bigger droplets 
being closer to the wall than the smaller droplets. Consequently, wall interactions may not be negligible, espe- 
cially for the larger droplets. The analysis by the Merritt and Subramanian (1989) shows that for an isolated 
bubble that is stationary under the action of buoyancy and thermocapillary forces, wall interaction effects 
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tend to move the bubble in the direction of the temperature gradient. That is, the presence of walls appears 
as though the thermocapillary forces are enhanced (or the buoyancy forces are reduced), up to separation dis- 
tances of even 10 to 15 radii. No theory seems to exist for the case of droplets. Our experiments seem to indi- 
cate that precisely such a phenomenon may be occurring; (iii) natural convection effects have no 
contributions to the force balance on the droplet in linear regimes as the temperature distribution in conduc- 
tive heat transfer in the vicinity of the droplet is symmetric about its mid plane. However, when convection is 
important, as in the present experiment, natural convection in the host fluid is expected to play a role. This 
is not included in the linear theory. Once again, since we have not measured the temperature distribution 
near the droplet we are unable to evaluate whether the contribution from natural convection enhances or 
reduces the migration of the droplet; and (iv) lastly, the effects of mass transfer for short times are unknown 
and may have an influence; mass transfer definitely influences the droplet motion for long times. 

In view of the lack of a theory which also considers the convective terms in the energy equation, an at- 
tempt has been made to describe a unified empirical fit to represent results of these experiments. Applying a 

1 dT 

multiple linear regression (Neville and Kennedy, 1964), using — and as two independent variables, a fit 

has been calculated for the three different temperature gradients and droplet diameters, and they are shown 
in Figure 7. The empirical equation is expressed as follows 
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Te *~ Tm - k [{dY] 


( 4 ) 



b 


where a = -0.7903 andb = 1.639. The linear theory curves for droplets of 3 and 9mm sizes are also 
plotted in the same coordinate system for a comparison. The value of the theory of Young . could be ap- 
proached with decreasing temperature gradient and increasing the droplet size. Again, the trend is qualita- 
tively the same, but different in magnitude. 

While attempts have been made to calculate the thermocapillary migration of bubbles in zero gravity con- 
ditions (Szymczyk etai, 1987, Shankar and Subramanian, 1988, Balasubramaniam and Lavery, 1989), there 
do not appear to be any studies under normal gravity for droplets for non-zero Reynolds and Marangoni 
numbers. These calculations would have to be carried out to compare with the experimental results. 

It is also well known that the motion of drops and bubbles is affected by surfactants which can adsorb on 
the interface and influence the interfacial tension coefficient (Harper, 1972). In experiments with water, 
Thompson et al . (1980) attribute the lack of motion of gas bubbles in a temperature gradient to the possible 
contamination of water by surface active impurities. The above consideration prompted a flow visualization 
experiment using droplets of water-methanol mixture (W-M, 50-50 by weight) in the 5 CS. silicone host-liq- 
uid. Using laser flow visualization, the thermocapillary flow inside the droplets were compared for W-M and 
vegetable oil, for comparable temperature gradients and droplet sizes. Even though the vegetable oil is 
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much more viscous than the water-methanol mixture, a flow pattern (Figure 6) resembling Hill’s spherical 


vortex (Harper, 1972) was easily seen in vegetable oil droplets. In contrast, both the surface flow and the 
return flow in W-M seemed to be considerably retarded. 

5 Conclusions 

Thermocapillary flow inside and on the interface of immiscible droplets in a host fluid have been, for the 
first time, visualized with directions consistent with the sign of the measured thermal interfacial tension coef- 
ficient (i.e., < 0). The flow persisted for long times. For short times after injection, droplets migrate 

toward the hot upper wall. Mass transfer affects the migration of the droplets for longer times. Our ex- 
perimental results follow the linear theory predictions qualitatively, but are quantitatively different. The 
droplets migrated farther away from the density-matched temperature compared to linear theory predic- 
tions. This appears to be not because of stronger thermocapillary flows, but due to reduced effects of buoyan- 
cy. We speculate that this may be due to (i) effects of convective heat transfer as the Marangoni numbers are 
non-negligible, (ii) wall interaction effects, (iii) natural convection effects which have not been accounted for 
in the linear theory, and finally (iv) the effects of mass transfer for short times that are unknown. Flow 
visualization studies in droplets of water-methanol mixture in silicone oil have revealed a significantly 
retarded flow pattern, compared to vegetable oil droplets. This suggests possible surfactant influence in the 
water-methanol mixture. 
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Based on our experience, reduced gravity experiments are perhaps the only way to provide direct and un- 


ambiguous migration data for non-negligible Marangoni numbers. 
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Table 1 

Properties of Liquids 


\Property 

p g /cc 

- 1 

V cs* 

cal 

1 

cm ~ 

Pr 

Fluid 

at P 

= 15 C 

cm sec °C 

sec 

Silicon 



0.00026 


71.75 

Oil 

0.911 0.0001197 

5 

0.0007 

Vegetable 

Oil 

0.91 1 0.000703 

66.9 

0.0003766 

0.00077 

868.83 


Table 2 

Measured Data and Comparison with the Linear Theory Predictions 


dT_ 

dz 

a 

_C 

cm 

D 

mm 

V m 

MEASURED 

f i m /scc 

V P 

PREDICTED 

LX 

p m / sec 

Vm 

Vm D 
2 a 

da d T n 2 

i > 7f-d 7° „ 

An 1 

23.63 

8.0 

143.8 

390.4 

2.72 

8.25 

1049.3 

19.81 

7.6 

105.4 

311.0 

2.95 

5.75 

793.9 

14.35 

7.63 

84.0 

226.2 

2.69 

4.60 

579.9 

23.63 

8.0 

143.8 

390.4 

2.72 

8.25 

1049.3 

23.63 

5.45 

125.6 

266.0 

2.12 

4.90 

487.0 

23.63 

4.1 

93.1 

200.1 

2.15 

2.74 

275.6 
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Experimental Apparatus 
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Figure 4 
Raw Data 
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Figure § 

Comparison of Present Results with Linear Theory Predictions 



Figure 6 

Long-Time Exposure Photo of Flow Inside a Droplet 
(D - 9mm, dT/dz - 6 9 C/cm) 
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Figure 7 

Unified Curve-Fitted Results and Comparison with Linear Theory 
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